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ABSTRACT 

o 

' We derive a relativistically covariant (although not manifestly so) equation for the 

distribution function of particles accelerated at shocks, which applies also to extremely 
relativistic shocks, and arbitrarily anisotropic particle distributions. The theory is for- 
mulated for arbitrary pitch angle scattering, and reduces to the well-known case for 
small angle scatterings via a Fokker-Planck approximation. The boundary conditions 
for the problem are completely reformulated introducing a physically motivated Green's 
function; the new formulation allows derivation of the particle spectrum both close and 
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{NJ ' far away from the injection energy in an exact way, while it can be shown to reduce to 
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£N) ■ a novel way. Contact is made with the Newtonian treatment. 



. Subject headings: shock waves - cosmic rays 

1. Introduction 

h 

The theory of particle acceleration at shocks is currently unsatisfactory. We cannot follow this 
process from first principles, i.e., the development of a collisionless shock, with the ensuing injection 
of non-thermal particles, nor do we understand the later shock evolution, which should include the 
interaction between the population of high-energy, non-thermal particles, and the shock structure. 
Even in the test-particle limit, to which this paper is totally confined, there are at least two major 
problems. First, the properties of the scattering agents are very poorly known; this problem will 
not be dealt with in this paper. Second, the general foundations of the theory have not been laid 
down. In particular, there is no general equation for the particle distribution function (DF) for 
arbitrary (i.e., even relativistic) shock speed, and no general treatment is available when the DF is 
anisotropic. Most results available to us either have been derived in the vanishing speed (thus fully 
isotropic) limit (Bell 1978, Blandford and Ostriker 1978), or descend from numerical (Bednarz and 
Ostrowski, 1998) or semi-numerical (Kirk and Schneider 1987) approaches. This is true even in the 
idealized situation of a plane shock expanding forever in a uniform medium, with particle scattering 
due only to pitch-angle scattering. It is the aim of this paper to derive a sufficiently general equation 
describing in a relativistically covariant way this process, to elucidate the conditions leading to the 
determination of the particle spectrum, and to apply these results to the extremely relativistic 
limit, where the DF anisotropy is expected to be maximal. In passing, a number of old results will 
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be recovered: it will be shown that the probability of reaching infinity is independent of the particle 
momentum, that the spectrum is the superposition of many bumps, each corresponding to a set of 
particles which have crossed the shock 0, 1, 2, N, .. times, and that, for very large particle energy, 
the superposition of these bumps spectrum leads to a power law in the particle impulse: all these 
results will be shown to hold for all shock speeds, Newtonian, relativistic or intermediate ones. 
Also, the explicit dependence upon the injection spectrum will be presented. In a future paper, 
this formalism will be applied to the extremely relativistic limit, including the effects of small or 
large angle scatterings, and a mean magnetic field. 

The plan of this paper is as follows. In Section 2, we derive the new equation for the particle 
DF, in an explicitly relativistic covariant way. In Section 3, we reformulate the boundary conditions 
to which this problem is subject, in a physically motivated way, finding eventually the dependence 
of the distribution function upon the injection spectrum. In Section 4, we concentrate on the 
distribution function at energies large compared with those of injection; here we show that it must 
be a power law, showing what fixes the power law index analytically, and deriving the relativistic 
generalization of Bell's law. We also make contact with the Newtonian limit. The important 
probability distributions P<j and P u , whose existence is merely postulated in the previous sections, 
are explictly written out in terms of eigenfunctions of the angular part of the scattering equation 
in Section 5. In Section 6, we summarize our results. 



We consider a shock propagating at arbitrary speed in an homogeneous fluid. We place our- 
selves in the shock reference frame, and choose coordinates such that z = identifies the shock 
position. The upstream fluid is located at z < 0, so that the fluid speed is always > 0, both up- 
stream and downstream. Suitable jump conditions at the shock are assumed (Landau and Lifshitz 
1987), but they will not be necessary in this section. 

We assume now that we can break down the effects of the magnetic field into two parts, as is 
customarily done. The first part is due to a small-coherence length magnetic field, possibly self- 
excited by the particles; this provides for an effective scattering, which is included in the collisional 
term. Another part, however, is due to a long-coherence length field, which provides instead for 
a smooth deflection of particles in phase space; this effect is included in the convective term (see 
below). We are thus using a dichotomic description of the magnetic field; most of the results to 
be described here, which are wholly unavoidable within this description, ultimately depend on the 
correctness of this assumption on the magnetic field. 

The distribution function is of course subject to the collisionless Vlasov equation, supplemented 
with suitable collisional and injection terms. In somewhat abstract form, we have 



2. An equation for the distribution function 
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The right hand side is the collisional term, with collisions per unit proper time, plus the rate of 
particles injected at the shock. 

The convective term is 

ds ds dx u ds dp u ' ^ ' 

In the above equation, the term containing the derivative with respect to the energy is not usually 
explicitly written out, imposing instead ab initio that particles be on their mass shell; here we 
do this instead (showing later on that it makes no difference) in order to emphasize the covariant 
character of the equation. In fact, both the first and the second term on the right hand side are 
invariant for Poincare group transformations. We make use of this by computing each term in 
different frames of reference, and with different orientations for the coordinate axes. In the first 
term, we place ourselves in the shock frame and orient the z axis along the shock normal, so that 
we can make use of d/dt = d/dx = d/dy = 0, because we assume the shock to be in steady state, 
and to have planar symmetry. We also write the z component of the particle four-speed in terms 
of quantities in the fluid frame, 

dz 

— = lpl {u + n) , (3) 

where [i is the cosine of the angle which the particle speed makes with the shock normal, z, in 
the fluid frame, u and 7 are the shock speed (in units of c) and Lorentz factor speed with respect 
to the fluid, and 7 P is the particle Lorentz factor also with respect to the fluid. All particles 
are here assumed relativistic in any reference frame (so that their three-speeds are always ss c), 
even those which are just injected. Strictly speaking, this is accurate only for relativistic shocks. 
The generalization to include Newtonian particles, though easy in principle to do, introduces some 
cumbersome formulae (with radicals) which complicate this work uselessly: it will not be considered 
any further. 

The second term on the right hand side of Eq. 2 will be instead computed in the fluid frame, so 
that we may neglect any motional electric field. In this case, the particle energy is conserved, and 
the term including the derivative df /dE disappears. The remaining components can be written as 
(Berezinsky et al, 1990, Ch.IX, Eq. 9.16) 

dp^di = Md^di _ d£ 

ds dp" ds dt dp" lpUJ d(j) 1 ' 

where u> = eB/E is the particle Larmor frequency in the fluid frame, and <f> is the longitudinal 
angle around the direction of the magnetic field, in momentum space. 

To estimate the collisional term, we first place ourselves in the fluid frame, and then remark 
that, in full generality, such term has the form 

(%) = ~ d ^ )f + j ^ ^ f ^' ■ (5) 

Here, we have made use of the assumption of linearity: we assume, that is, that the magnetic field 
turbulence responsible for the scattering is independent of the particle density. We consider particle 
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scattering when there may be a correlation between the directions before and after the scattering 
event, so that the scattering coefficient w(fi,fi') depends on both. Also, we are not assuming 
the diffusion approximation: the above equation applies to arbitrarily few, and arbitrarily large- 
angle, scattering laws. It would reduce to the well-known Ginzburg-Syrovatskii (1964) form in the 
diffusive approximation. The first term on the right hand side represents scattering of particles 
away from their direction of motion, while the second term represents the total contribution to the 
distribution function in the direction of motion under consideration, by scattering from all other 
directions. Clearly, d and w are related through 

d{n) = J ^ w(n',fj,)4i' (6) 

which simply expresses probability conservation, w (and thus also d) may depend upon p. Often, 
one makes the simplifying assumption that the p-dependence in w (and thus in d) factors out; we 
shall comment later on why this may come in handy, but will not assume this. We also remark 
that w, d > 0. 

Since we need df/ds, while above (Eq. 5) we have evaluated the collisional term in the fluid 
frame, we correct it by means of dt/ds = j p . Also, we wish to include the injection term, 4>i n j. 
We assume that injection takes place at the shock only, but we leave the energy and direction 
dependence of the term otherwise arbitrary: (pi n j = Gi n j(p, /i)5(z) where 5(x) is Dirac's delta. 

Putting together all of the above, we find 

7(« + = + J wQjl^W)^') +u^ + G inj (p,(i)8(z) , (7) 

which is the equation we were searching for. It is worth remarking that all terms in the above 
equation are computed in the fluid frame, except for the space dependence, and its attaining 
derivative. Also, the above equation is valid both upstream and downstream, but now care must 
be taken because the quantities p, fj,, fi' ... are all evaluated in the local fluid frame, and thus 
quantities bearing the same name in the two distinct reference frames are not identical: they are 
related instead by a Lorentz transformation, which we give here for future reference. Using subscript 
u and d to denote quantities in upstream and downstream frame respectively, he have 

Pu = Pd7r(l - u r fi d ) ; fi u = — — , (8) 

1 - U r fi d 

where u r and j r are the modulus of the relative speed between the upstream and downstream 
sections, and its associated Lorentz factor. 

Whenever the scattering in each individual event is small, it is possible to apply to the above 
equation the very same treatment which leads to a Fokker-Planck equation. In this case one obtains 
(Kirk and Schneider 1987) 

1{U + ^J- Z = c| "^ 2) §) + G inj (p^)6(z) , (9) 
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where the effective diffusion coefficient is related to the r.m.s. angle of diffusion. It is worth 
remarking that it may occur that the above equation holds for the downstream section, but it is 
not justified for the upstream section, where the full Eq. 7 has to be employed. The reason is that 
the diffusive approximation on which the above equation is based requires each scattering event to 
lead to a r.m.s. deflection 0(1/7), where 7 is the fluid Lorentz factor with respect to the shock. 
When the shock becomes extremely relativistic, thus, the correctness of the diffusive approximation 
breaks down: it follows that one may use Eq. 9 for the downstream section, but is forced to use 
Eq. 7 for the upstream one. 

In the following, we drop the term due to the long-coherence length magnetic field, which will 
be considered in Paper II. 



We begin by remarking that in this problem there is a net mean flux of particles across any 
surface; if the surface is steady in the shock frame, then the flux is also time-independent. Since 
the particles' speed is always assumed « c (« 1 in our units), the infinitesimal flux across a unit 
shock area, per unit time, is given by 



depending upon whether we choose to express / in terms of variables in the shock, upstream, 
or downstream frames. The representation in terms of downstream variables is especially useful, 
because in the downstream section the downstream momentum pj is conserved, so that a net flux 
across the shock of particles of impulse pd can be due only to the fact that some fraction of these 
particles is advected to downstrean infinity, with probability given by Eq. 21. 

From now on, in this section, we shall use coordinates in the downstream frame exclusively, 
which we thus indicate, for ease of notation, without subscripts. 

The mathematical treatment of this problem must be physically motivated: we envision a 
system where particles are injected at low energies at the shock, from which they cannot reach 
upstream infinity, because they would have to swim upstream, but from which they do reach 
downstream infinity. In order to study how particles generated at the shock reach downstream 
infinity, we first concentrate on the downstream part of the problem. 

Consider now a downstream section of finite length L, beginning at the shock; then place two 
observers at either end, and have them inject particles into the downstream section. Then they 
may collect outgoing particles at either end. The problem is obviously physically well defined. The 
speed of a particle (in downstream frame coordinates) with respect to the shock is 



3. Obtaining the particle spectrum 



dJ = HsvlfdHsdPs = lu{u + ^u)vif<^udpu = ld{ud + ^Pdf&ddPd 
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Thus, the observer located at the shock injects particles for fi > — Ud (here Ud is the downstream 
fluid speed with respect to the shock), and collects the outgoing ones (i.e., those with \i < —Ud). 
So reasonable boundary conditions at the shock can be posed for fi > —Ud only. The observer at 
the other end of the downstream section disappears as L — > oo, so that there can be no entering 
flux there; we can only impose the regularity condition that / does not grow to infinity, nor goes 
to zero. 

Similarly, for the upstream section we see that we can impose boundary conditions for fi < —Ud 
at the shock; at upstream infinity, the only reasonable boundary condition is that /—►(), because 
we cannot expect our particles to swim against the fluid advection all the way to upstream infinity. 

We now see that the boundary conditions for the upstream problem are provided by the out- 
going particles of the downstream section, and vice versa the outgoing particles of the upstream 
section provide the boundary conditions of the downstream section, thusly setting up an obvious 
(fixed point) problem. In order to make this explicit, let us introduce the conditional probability 
Pd{^ini Hout)dl JL out that a particle, given that it crossed the shock toward downstream along a direc- 
tion Hi n , will recross it along the direction fi ut- Then, calling Jj n = (ud + Hin) f (l^in) the entering 
flux, the outgoing one will be given by 



)Jin(Hin) + G(p, l^out) = Pd* Jin + G ou t (12) 

J-u d 

where we have also included in the outgoing flux the part deriving from the injection G ou t = 
Gi n j(nout) for n < —Ud- The symbolic notation will be helpful in the following. 

Before specifying the equivalent relationship for the upstream part of the fluid, we pause to 
establish the connection between Pd and the solutions of Eq. 7 (or Eq. 9) . Given the definition of 
Pd, it seems obvious enough that Pd is the particle flux (u + n)f(n) given by the solution of Eq. 
7 (or of Eq. 9) with boundary condition (ud + fJ-)f = S(fi — fj,i n ) (unit entering flux) for fi > — Ud 
at the shock (z = 0), and / regular at downstream infinity. In fact, we can see that this solution 
represents a monochromatic (in energy) flux of particles, all moving initially in the same direction, 
subject to pitch angle scattering, which is in turn responsible for kicking them back toward the 
shock, occasionally, along a direction fi ut, with nonuniform probability. 

One way of looking at Pd is by realizing that it is somewhat related to the Green's function 
for a diffusive problem. As an analogy, consider the paradigmatic diffusion equation 

I - <»> 

which is solved by convolving the given boundary conditions with its Green's function 

2(vrz)- 1 / 2 cxp(-(/x - ^f/Az) -> 5(ji - (14) 



as z — ► 0. Here too, the Dirac's delta appears in the boundary conditions, and is smoothed 
out as the time variable (z) evolves. The major difference is that suitable boundary conditions for 



- 7- 



the above equation include the whole range in fi, at z = 0, while in our problem we can specify 
boundary conditions only for fi > — u^- The function P^ is the projection of the solution of the 
problem on exactly that part of the boundary at z = where we cannot specify the boundary 
conditions, a part which does not exist for Eq. 13. Still, because of this relationship, we will refer 
to Pd and to P u , shortly to be determined, as the Green's functions of the problem. 

As an aside, we show that, under some physically realistic conditions, P& does not depend upon 
the particle impulse p. This may be at first surprising because both coefficients d(/u,) and w(p, //) 
in Eq. 7 are allowed to depend on p: the proof follows. Rewrite Eq. 7 outside the shock (which 
means that we can drop the injection term) without the term depending upon the magnetic field, 
using as a new variable y = d(fi)z; this can always be done because, as noticed above, d > 0. We 
obtain 

7 (u + /x)^ = -/ + J g(v^')f(p'W . (15) 

Since we assumed that the p-dependence in w and d factors out, we see from Eq. 6 that g is 
independent of p. Then from Eq. 15 we see that p has disappeared from the problem altogether, 
and that any solution can depend on p only through the term y = d(jx)z. However, Pa is the flux 
at z = 0, the shock position, so that we see that when this is computed (y = 0), all dependence 
on p drops out: P& does not depend on p. This result was known to Bell (1978), who proved it for 
Newtonian shocks; the proof given here is however valid for all shock speeds. 

Upstream, we shall also have a probability distribution P u (Houti fJ-in) that a particle leaving 
downstream along a direction p ou t will recross the shock along a direction p,i n . An argument entirely 
identical to the one given above shows that, often, P u does not depend upon the particle's impulse. 
There is however a small complication, in that particles of impulse pi in the downstream frame, 
exiting it along a direction p out and re-entering it along a direction fii n , emerge with a different 
impulse, p (as measured again in the downstream frame) given by 

1 — U r H ou t 

V = -t~Z Pi = Gpi ■ 

G is the energy amplification a particle receives, as it cycles once around the shock. Since — 1 < 

H-out < ~Ud, and -u d < fi in < 1, we have 

1<G<\±^. (17) 
1 — u r 

Thus, particles leaving the downstream section with impulse pi provide the reentering flux at a 
different impulse. A trivial computation yields the entering flux at impulse p as 

JiniPi fJ'in) = / dp ou tPu{pouti A*in) f^j ' ) Jout{Pi: l^out) 

J-l \l-U r fl out J 

+G in (p, p in ) = P u -k j' out + G in (p, p in ) (18) 

where again we added to the incoming flux the injected one Gi n = Gi n j(^i n ) for pi n > —u^, and 
the prime over J out reminds us that the outgoing flux is to be computed at particle energy pi given 
byEq. 16. 
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A comment on normalization is in order. All particles that enter upstream will in due time be 
brought backwards across the shock by the fluid's advection, so that we expect 

/ dfJ, in P u (fJ, out , ftn) = 1 (19) 
J -u d 

while the analogous statement for the downstream Green's function P^ does not hold for the same 
reason: a fraction of all particles will be advected toward downstream infinity, and the returning 
flux will be smaller than the departing one. In general, thus, f_^ d dftoutPdil^in, fJ-out) < 1- The 
average probability of coming back to cross the shock is given by averaging Pj over the whole 
incoming flux: 

ini l^out) Jin\^in) 
J~Ud d/J'inJin([J'in) 

When the contribution of injection can be neglected, we see from the above, and from Eq. 12 that 



/_ \ d d/J-out Joutifl' 



out I 



P = . (21) 

in) 

The reason for this is quite simple: in a steady state, like that assumed here, there can be no 
accumulation of particles except at infinity; since there is a larger flux entering the downstream 
region than leaving it, there must be an accumulation of particles in it, which can thus be realized 
only at downstream infinity, i.e., by leaving the region of the shock. 

Our cardinal equations, Eqs. 12,18, can now be simply combined to obtain 

Jin = Pu* Pd* JL + (Pu * G' out + Gi n ) = Q * j' in + X . (22) 

This equation can be solved iteratively, using the term in parentheses as the first guess: 

Jin,(0) = X i Jin,(l) = Q * X> + X , Jin,(N) = Q * J'in,(N-l) + X ■■■ (23) 

to arrive at the solution in the form 

J in = X + Q*X' + Q*Q*X" + Q*Q*Q*X'" + .... (24) 

The interpretation of this equation is simple: X is the flux of injected particles, including both 
those that are directly moving in and those that were initially moving out, but have been turned 
back exactly once, The operator Q transforms an ingoing flux into another ingoing flux, made 
of the fraction of all particles that have completed exactly one loop around the shock. Thus the 
terms Q * Q * X" and so on represent the particles which have made 2, 3, ...N, ... loops around the 
shock, after injection. The fact that this theory manages to recover the well-known fact that the 
spectrum is the superposition of infinite bumps, each corresponding to particles that have looped 
an integer number of times around the shock, can be considered a useful test which has been 
successfully cleared. It also has the advantage of yielding the distribution function at the shock for 
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the first time both very close to the injection energy, and very far from it. Also, it substantiates 
physical intuition that the solution of the problem must depend only on the injected flux (X) and 
on the scattering properties of the medium. Lastly, the equation concerning J ou t carries no new 
information: once Jj„ has been found, it can be plugged into Eq. 12 which yields a J out which 
automatically satisfies Eq. 18. 



4. The large-p limit 

We see from Eq. 16 that after half a cycle (only the upstream section changes the particles' 
energies), the energy of a particle just injected with energy pi n j may span the range (— 1 < 
mUout < -u d < n in < 1) 

1 + u r . . 

Pinj < P < Pinj~^ ~ = Prnax ■ \*"J 

This shows that the evolution of that part of the particle distribution function which describes 
particles making many loops around the shock, must extend toward larger and larger energies; in 
fact, since the minimum amplification is G = 1, the process under discussion does not describe a 
spreading of the injection distribution around the initial value, but a drift towards larger and larger 
energies (plus a spreading, of course): Fermi acceleration of Type I. So the question naturally arises 
of what is the asymptotic spectrum of particles, as we consider asymptotically larger energies. 

If injection provides particles of greatest energy Pi n j, we may neglect the injection terms in the 
above equations for energies exceeding p ma x, and we thus consider Eq. 22 in this case: 

Jin(p,fj>)=[ df-outPuifj-out, n) ( ~r — f (%P d (£, fj, out ) J in {\ — Ur ^ m p, g) . (26) 

J-l V-L - UrfJ-out J J-u d L — UrHout 

We show now that this equation, except for the null solution, has a solution which is a simple 
power law in p. Let us first consider two situations, where injection occurs at slightly different 
energies, po and po + 5p$. At large energies, p > p max , the spectra of these two situations will 
differ by very little, d,Ji n /dpo 6po. The difference between two solutions of the inhomogeneous 
problem (i.e., including injection) must be a solution of the associated homogeneous problem (the 
one without injection, Eq. 26), so that dJi n /dpo is a solution of Eq. 26. For obvious dimensional 
reasons, Jj„ = h(p/po, [i)/p 3 so that (using h(z,fi) = dh/dz) 

dJi n 3h h 

+ -5— (27) 



dp p 4 p 3 po 



dJin 



and, upon eliminating h, 



dp p 2 Pq 

:--3J -1 
op PO Opo 



(28) 



Jijhn = _ 3Jm _ ]_dJn (29) 
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which shows pdJi n /dp to be the linear combination of two solutions of the homogeneous equation 
26 (for p > Pmax), and thus a solution itself. If we now assume Eq. 26 to have a unique solution, 
we find that pdJi n /dp must be proportional to Jj n , when injection can be neglected. We thus have 

' (Um (p,ri = -- (30) 



Jin(p,V) dp ' p 

where s is a constant yet to be determined, and the obvious solution 



J in (r,ri= iud+ £ )M (3D 

follows, where we have arbitrarily factored out Ud + \i to remind us that J is a particle flux; thus, 
g{pL)/p s is the particle distribution function. It must be remarked that this power-law dependence 
has been derived without assuming that and P u are independent of the particles' momentum p. 

The importance of the above is obvious: it shows that the particle distribution function is 
a power law in the particles' momentum, for every shock speed. This result had been obtained 
before in all numerical solutions for arbitrary shock speeds, and in analytic form in the Newtonian 
limit, but in the full relativistic regime it was put in by hand (Kirk and Schneider 1987), and 
the boundary conditions could have equally well been satisfied by any other function with a free 
paramemter. The present derivation is universal, and it has the rather unexpected result (at least, 
for this author) that the power law solution is exact, not asymptotic, as soon as we consider energies 
exceeding the injection ones, i.e., for p > p max , and not just p 3> p ma x- 

With the result provided by Eq. 31, our principal equations, Eqs. 12 and 18, assume the 
following form, when injection may be neglected: 



(u d + f i)g(ri= f dUPd&riiua + OgiO (32) 
(u d + n)g{n) = / d(i out P u (fj, out , //) — (u d + Hout)g(^out) ■ (33) 

J -l \l-Ur^outJ 

The first one concerns the outgoing flux (with respect to the downstream section) , and thus is valid 
for — 1 < (i < —Ud, while the second one concerns the entering flux, and applies for — Ud < u < 1. 
Also, use of Eq. 31 in Eq. 21 shows that the probability of reaching downstream infinity does not 
depend upon the particle's momentum, p. 

Again, they can be combined to obtain 

(u d + fi)g(u)= f 1 dtQ T (Z,fi)(ud + Og(0 (34) 

J— Ud 



u d /i .. ,.\ 3-s 



where 

Q T (tu) = , 

l_i \ 1 — u r v 



Q r {M= I <bP u {v^)Pd{^)[^-^) . (35) 
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which applies to the incoming flux only, — u d < fJ, < 1. Once again, the outgoing flux carries no 
new information: if the above is solved, plugging the solution into Eq. 32 returns an outgoing flux 
which automatically solves Eq. 33. All of the problem's information is in one of the subintervals 
only, even in the case of the homogeneous equation. 

The above equation, for an arbitrary value of s, does not have a solution. The equation which 
does have a solution is 

X(u d + ^g( f i)= f 1 c£Q T (£,/x)(u d + 0<7(0- (36) 

J -u d 

In fact, the above equation is an homogeneous equation of the Fredholm type, with smooth and 
bounded kernel over a compact interval (Courant and Hilbert 1953), which has at most a finite 
number of solutions, each belonging to an eigenvalue A which is to be determined simultaneously 
with the eigenvector g(fi). We thus see that what fixes s is the requirement that A = 1: the physical 
value of s is that which has a unitary eigenvalue. We are forced by our limited mathematical skills 
to assume that Eq. 34 has exactly one solution, neither more nor fewer, but, apart from this, the 
task of determining the particle distribution function is completed: Eq. 36 with A = 1 yields the 
particle distribution function at the shock, and the all-important particle spectral index s, even in 
the case of arbitrary particle anisotropy. 



4.1. The Newtonian limit 

In order to gain physical insight into the condition that fixes s (A = 1), we rewrite Eq. 36 
somewhat. First, we write Eq. 33 as 

X(u d + n)g{n) = I (faoutPuiVout, A*) i : — {ud + Hout)g(Hout) ■ (37) 

J-l V-L - UrUoutJ 

This can now be integrated over the whole range in fi, —Ud < fi < 1, and divided by the whole 
outgoing flux, to obtain 

1 — ( 1— \ 3— s 

j\ d dfi(u d + mMaO }- Ud <bLi d 4*outPu(pout,v) ) (u d + AWMAW) 

A- 



J_" d 4l OU t(ud + Hout)g(Vout) f_i d d^out{u d + Hout)g{Vout) 

(38) 

The term on the left hand side can be rewritten by means of Eq. 21 as X/P. The term on the right 
hand side is clearly the average value of the (s — 3)-th power of the amplification G (Eq. 16), over 
the whole particle distribution: < G s ~ 3 >. The above is thus 

X = P < G s ~ 3 > , (39) 

and demanding that A = 1 implies that 

P < G s - 3 >= 1 . (40) 
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In the Newtonian limit, in which G — 1 <C 1, we obviously have < G s 3 >^< G > s 3 , and thus 
the energy spectral index k = s — 2 is given by 

logP 

s - 2 = k = 1 - - °% (41) 
log < G > 

which is exactly Bell's equation for the spectral index. We have thus established that the condition 
A = 1 is equivalent to P < G s ~ 3 >= 1, which appears to be the correct relativistic generalization 
of Bell's Newtonian results. 

In the Newtonian limit, where g{ft) ~ constant, independent of [i (the isotropic limit), it is 
easy to derive from Eq. 21 that P ~ 1 — u, and from Eq. 38 that < G >~ 1 + u, from which s = 4, 
as per Bell's results for a strong shock. 

There are some differences between the fully relativistic approach and the Newtonian one. 
First, in the relativistic approach in general < G s ~ 3 > cannot be approximated by < G > s ~ 3 , 
contrary to Peacock's (1981) arguments. Second, here G and P are functions of s, so that Eq. 
40 sets up a transcendental equation for s. Third, Eq. 40 is only of symbolic value, since it 
requires knowledge of the anisotropic particle distribution function, which is not known a priori. 
But luckily, the full relativistic approach developed here returns simultaneously the anisotropic 
distribution function (Eq. 34). 



4.2. Bits and pieces 

A useful identity is obtained by rewriting Eq. 18 without the injection term as 

Jin{Piftin) f d , D / n Jout{Pii ft out) fAO^ 
<jftout" u \ftout, ftinJJZ To , 



- U r m n ) 6 J_ x (1 - UrUout) 6 

integrating over using Eq. 19, and then inserting Eq. 31, to obtain 

/i r-ud 
dft in (l - U r fl in ) S ~ 3 (u d + ftin)g(ftin) = / ~ «r ftoutY' 3 {u d + ftout)g(ftout) 

-u d J-l 

which of course is not trivial, because the intervals of integration are different. 



(43) 



It is worth noticing that, in all of Section 3, no reference whatsoever has been made to the 
explicit form of the right hand side of Eq. 7; all that is necessary is to determine the equation's 
Green's functions P& and P u , but otherwise all results apply to all right hand side members, and 
can then be seen as a property of the left hand side member (u + ft)df jdz. They are thus of very 
general validity. 

Lastly, we have assumed that Eqs. 34 and 36 has one, and just one solution. If it had none, 
then the steady state problem would have no solution, and the problem would be intrinsically time- 
dependent. If it had several, then a rather interesting case would arise, whereby shocks moving 
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at the same speed in the same media, but with different injection properties would have different 
asymptotic spectral indices s. The near universality of the index s points of course in the opposite 
direction, but it is also true that a proof of the existence of a unique solution is beyond our modest 
abilities. 

5. P u and P d 

In this subsection we show how to build P u and P d . For sake of definiteness, we concentrate 
on P u , but a wholly analogous treatment holds for P d . 

We now wish to find an explicit expression for P u (fi ut, Hin)- That this quantity exists, it is 
physically obvious: it represents the probability that a net flux of particles leaving the downstream 
section along a direction fj, ou t (in downstream variables), reenters along a direction fj,i n . We can 
imagine an experiment to determine this quantity: an experimenter located next to the shock 
with a small cannon shooting out particles along the direction fi ou t, may collect them on a screen, 
determining how many come back along //j n . It corresponds to a solution of Eq. 7 or 9 without the 
injection term, with boundary condition 

dJ = 6(nd- Vout)du d = 6(n u - ^ t )da u . (44) 

(u) 

Here, n J t is the value of \x ou t in the upstream frame variables, and the last identity is a trivial 
property of Dirac's deltas. 

In order to find P u , we now use separation of variables to solve Eq. 7 or Eq. 9. For instance, 
writing f([i u , z) = A(z)B(fi u ), we find the solutions of 

1 dA . 1 d fn, wi 2x dB \ 

Calling B n the angular eigenvectors, we find a generic solution of the above as 

/ = ^2 a n ex.p(X n z/'y)B n (n u ) (46) 

n 

where the a n 's are coefficients which we must choose so as to satisfy the initial conditions, Eq. 44. 
We must thus have, at z = 0, 

7(« + A*u)/ = Yl a ^&)i u + VuhB n (fi u ) = 5{fi u - . (47) 

n 

(u) 

Obviously, the coefficients a n 's will depend upon the point where the delta is located, fi Q J t . 

Here there is a tricky, but important point. The above eigenvector problem is well-known to 
have solutions belonging to both positive and negative eigenvalues (plus the A = case). Eigen- 
vectors belonging to negative values of X n , of which there is an infinite number, are physically 
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ill-behaved at upstream infinity, as can be seen from Eq. 46, so that we surely have to restrict 
ourselves, in all sums above to the well-behaved eigenvectors, i.e., those with A ra > 0. But does 
the above equation then hold? The answer is yes, as we now show. 

The fact that an infinite set of functions can be arranged to satisfy the above equation is called 
by mathematicians completeness, indicating that any function can be written as the superposition 
of BnS with suitable coefficients. If the above equation holds, in fact, multiplying both sides by 
any function F(fi^ t ) and integrating over the whole range of n^ t , we find that 

F ^u) = E (/ An(^ t )F(^ 7 (« + Hu)B n (n u ) , (48) 

proving that we can write any function as the superposition of functions from the set, with suitable 
coefficients, and thus that the set is, by definition, complete. Now, the set of all eigenfunctions B n 
is well-known to be complete, but here is the rub: the full set includes both the eigenvectors which 
have A > (which thus are physically well behaved at upstream infinity) and those which have 
A < 0, which diverge at upstream infinity, and which we thus cannot use. 

The half set of the well-behaved eigenvectors is not complete over the whole range in fx, but 
here we are helped by an unusual property of the eigenvectors of this equation: it can be shown in 
fact (Freiling, Yurko and Vietri 2002) that the half-set of all well-behaved eigenvectors (i.e., those 
with A > 0) is complete in the restricted range — 1 < [i u < —u. This property (sometimes called 
half-range half-completeness) implies that, for fi u < —u, we can always write 

l{u + lJ- u )f = E a "( n + VuhB n (n u ) = 8{nu - /x£) (49) 

n 

where however now the prime over the sum reminds us that the summation is extended only to 
well-behaved eigenvectors. For fi u > —u, the above equation is not verified, and the sum will give 
some function, which we can now identify with our required P u . Thus, P u is the continuation to 
the rest of the interval (fi u > —u) of the sum above, which inside the restricted range fi < —u 
satisifes the above equality. 

One word of caution should be stated here: the property of half-range half-completeness has 
been proved for the eigenvectors of Eq. 9, but not for those of Eq. 7. Still, since they are known to 
be so similar (both are called of the third, or polar type, one within the realm of integral equations, 
the other one of differential equations, to emphasize several formal similarities), and since the above 
physical argument makes it extremely plausible that P u exists also for this case, we shall make no 
further distinction between the two cases, and proceed as if this argument went through for both 
Eq. 9 and Eq. 7. 

We now find the coefficients a n 's. The well-behaved eigenvectors are not orthormal over the 
restricted range \x u < —u. They can easily made so by means of a well-known procedure (Schmidt's 
diagonalization, Courant and Hilbert 1953), which leaves the lowest order eigenvector unaffected 
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(except for the normalization), and suitably modifies the others in a finite number of steps. Once 
this new basis P n has been found, since it is complete (from the half-range half completeness), we 
easily find, for fi u < — u, 

(U + fX u ) £ Pn{&)Pn^u) = S(fl u - M2) (50) 
n 

and thus, by working backwards to Eq. 10, (for fi out < —Ud, fJ-in > —Ud) 

Pui^out, VinWin = {U + fi u ) ^ P n {^ t )P n {^ u )dlX u . (51) 

n 



6. Summary 

The major results presented in this paper are the following: 

• We have presented a relativistically covariant equation (Eq. 7) for the particle distribution 
function, which includes pitch-angle scattering (without assuming that the diffusive approx- 
imation holds), the effect of smooth magnetic fields, and particle injection at the shock; this 
equation applies separately in the upstream and downstream frames. 

• We have described the transport properties of the upstream and downstream media by means 
of two Green's functions, Pd and P u , which are independent of the particle distribution 
functions, and which have been explicitly constructed in terms of the eigenfunctions of the 
angular part of the problem (Eq. 51). 

• We have established two relations (Eq. 12 and 18) which provide the true boundary conditions 
of the problem, giving the flux entering the downstream section in terms of the one leaving 
the upstream section, and vice cersa, plus the injection terms. This naturally sets up an 
equation for the particle flux, Eq. 22. 

• We have solved this equation, Eq. 24, in a general way, which yields the particle spectrum 
both close to and far from the injection energies. 

• We have shown that under physically realistic conditions, the probabilities of returning to the 
shock, P u and Pd, do not depend upon the particles' energy p. 

• In the limit of large particle energies, we have shown that the spectrum is a pure power law 
in the particles' momentum, even when the probabilities of returning at the shock, P u and 
Pd, depend on p (Eq. 31). 

• In this limit, we have then simplified our main equations, arriving at a system of equations 
(Eq. 32, 33, or their combination, Eq. 34), and further showing that the requirement that 
Eq. 34 has the eigenvalue A = 1 fixes the all-important asymptotic particle spectral index s. 
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• We have recovered the low speed, Newtonian limit of Bell (1978) and Blandford and Ostriker 
(1978), and provided a relativistic generalization (Eq. 40) which is however of little use 
because it requires knowledge of the particle anisotropic distribution function, which can be 
obtained only through the full solution of Eq. 34, which automatically incorporates Eq. 40. 

Though the treatment presented in this apper may appear rather abstract, it is not without 
practical consequences. For instance, the fact that the asymptotic spectral index s is not asymp- 
totic, but it applies as soon as p > p ma x = PinjiX + M r)/(1 — u r ), where u r is the modulus of the 
relative speed between the upstream and downstream fluids, implies that numerical simulations 
trying to compute s need not extend to very large particle energies, but can save precious com- 
putational time by sticking to energies exceeding p m ax 

by a small factor, dictated exclusively by 
numerical questions. Also, the treatment presented in Section 3 may be somewhat simplified by a 
judicious use of Fredholm's formulae, resulting in a considerable simplification for the calculational 
side of the theory. And lastly, they allow the treatment of the hyperrelativistic limit, as will be 
shown in a future paper. 
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